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We present an analytical theory on the heights of the quantum steps observed in the hysteresis 
loops of molecular magnets. By considering the dipolar interaction between molecular spins, our 
theory successfully yields the step heights measured in experiments, and reveals a scaling law for 
the dependence of the heights on the sweeping rates hidden in the experiment data on Peg and Mn4. 
With this theory, we show how to accurately determine the tunnel splitting of a single molecular 
spin from the step heights. 
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Crystals of molecular magnets, such as Fes and Mni2, 
have attracted much attention for their connection to 
macroscopic quantum tunneling and Berry phases [Q, H, |^, 
^. They may also have important applications in mag- 
netic memory and quantum computing jj, |^] . The earli- 
est and most spectacular observation on such a system is 
the quantum steps in the hysteresis loop of magnetization 
at low temperatures . 

These quantum steps are a manifestation of macro- 
scopic quantum tunneling, resulting from the tunneling 
between different spin states of large molecular spins 
{S = 10 for both Fcg and Mnia, S = 9/2 for Mn4). As in 
any many-spin system, this tunneling phenomenon is nat- 
urally complicated by the interaction between spins and 
other environmental effects. A full understanding of this 
phenomenon demands a successful theory on the quan- 
tum steps, describing their widths, shapes, and heights. 
So far, besides speculations, no serious attempt has been 
made to develop quantitative theories on any of the step 
features. 

In this letter we present a successful theory on the 
step height as the first step towards a complete theory 
on the quantum steps in the hysteresis loop. Since the 
height measures the tunneling probability between dif- 
ferent spin states, it is the most prominent feature of 
the quantum step, and holds the key to understanding 
the underlying tunneling dynamics. In Fig.|l|, we have 
adapted the experimental data on Fcg from Ref.[^, and 
show how the step height between spin states Sz = ±10 
changes with the sweeping rates. The data are compared 
to the Landau-Zener (LZ) model[0, which has been used 
to extract the tunnel splitting A of a single molecular 
spin from the step heights|^, ||. When we fit the LZ 
model with the data at the fast sweeping regime there 
is a dramatic difference at slow sweepings: a two-third 
suppression. 

By taking into account the dipolar interaction between 
molecular spins, our theory successfully gives the step 
heights measured in the experiment, as shown in Figj^. 
Two physical mechanisms influencing step heights are 
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FIG. 1: Comparison between the experimental data on FegQ, 
the LZ model, and our theory. The normalized step height 
5M/2Ms is the final fraction of up-spins Ffin after one sweep. 



identified: the spin shuffling in the evolving distribu- 
tion of dipolar fields and the dipolar interaction between 
spins in the resonance window. Furthermore, our theory 
reveals an a/A^ scaling law for the dependence of the 
heights on the sweeping rate a of the external field. This 
law is confirmed by the collapse of the experimental data 
in terms of the scaled sweeping rates (see Fig. As 
a direct application of our theory, we show how to accu- 
rately extract the tunnel splitting A of a single molecular 
spin from the step heights and the sample geometry. In 
our theory there are no adjustable parameters. 

We argue that nuclear spins do not appreciably affect 
the tunneling dynamics under a fast sweeping field except 
for a renormalization of the tunnel splitting. In the re- 
laxation experiments |^ where the external field remains 
constant, the tunneling is strongly affected by the nu- 
clear spins as recognized by Prokof'ev and Stamp [ p^ . 
However, in the sweeping field experiments, the role of 
nuclear spins is marginalized by the sweeping fields. This 
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can be clearly seen in Fig. 6ofRef. , where the relax- 
ation with a constant external field is shown to be much 
slower than the one with a sweeping field of the slowest 
rate appHed 0.04mT/s. 

We consider spin lattices, such as crystals of Fes, 
Mni2, and Mn4, in which the spins interact with each 
other through the dipolar potential, d(r) = E£){1 — 
3 cos^ 9)flQ / , where r is the displacement vector be- 
tween the spins, 9 is the angle between f and the 
easy axis, JIq is the unit cell volume, and Eu = 
'^{Sg^isY 1^0 gives the interaction strength. Our the- 
ory will be compared to the experiment, mainly on crys- 
tals of Fes where the experimental data on step heights 
are the most abundant . For simplicity, we focus on 
one step, that is, the tunneling between two spin states 
(for example, Sz = ±10 for Fes); it is rather straightfor- 
ward to extend our theory to study multi-step tunneling. 

We now have a system of Ising spins sitting at each 
site of a lattice. In a sweeping magnetic field along the 
easy axis, the spins will flip from one state to the other 
back and forth as a result of the tunneling driven by 
the sweeping. However, at any given moment, only a 
small fraction of spins are flipping by being in the reso- 
nance window while the others remain static. This can 
be understood by first considering an isolated spin in a 
sweeping field, which can be described exactly with the 
LZ model. In the LZ model the flipping occurs mainly 
in a tunneling time interval, when the Zeeman energy 
bias 7 = 2(7/iB5^o-ff@ between the two spin states 
caused by the changing external field becomes very small, 
I7I ^ A^i„/2. This tunneling time defines the resonance 
window, whose width Auji„ is the tunnel splitting A at 
the adiabatic limit and ^/2a {a = fi^) in the sudden 
limit p^ . Similarly for a spin interacting with other 
spins in a lattice, its resonance window is defined by 
I7 + Ci| ^ A.u,in/2 for spin i, where is the Zeeman 
energy of spin i caused by the dipolar field from other 
spins. Since the dipolar field felt by spins is a distribu- 
tion, only a small fraction of spins are in the resonance 
window at any given moment. 

With this physical picture in mind, we can write down 
the evolution equation for the fraction F of up-spins. If 
spins in the resonant window fiip with probability Pwim 
we have 



dF_ 



(l-2F)i?(-7,F)P^ 



(1) 



where D{^, F) is the normalized distribution of dipolar 
field with the fraction F of up-spins randomly located 
throughout the lattice, and j — at represents the sweep- 
ing field. The combination (1 — 2F)I?(— 7, F) is the differ- 
ence between the fractions of up-spins and down-spins in 
the resonance window. We want to solve Eq. (|l|) with the 
initial condition F = 0, that is, all the spins point down- 
wards at the beginning. The result i^/m = F{t — > 00) is 
the fraction of up-spins at the end of the sweep, or the 



normalized height 5M/2Ms of the quantum step between 
the two spin states. However, we need to first find what 
Pwin is, and how to calculate the distribution function 
Di^,F). 
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FIG. 2: The flipping probability obtained with the nonlinear 
LZ model. J/A = 0.0 corresponds to the hnear LZ model. 
The flipping probability is suppressed for positive J, and en- 
hanced for negative J. 

Without the dipolar interaction, the flipping proba- 
bility Pwin would be given by the LZ model, that is, 
Pyjin — — \ ^ exp(— ^^). With the dipolar interac- 
tion, it is no longer trivial to calculate Pwm- During the 
tunneling time defined by Au,i„/a, a spin inside the res- 
onance window feels two kinds of dipolar fields: one from 
spins outside the window, the other from spins inside the 
window and trying to flip together. The former remains 
static during the short-time flipping process; it merely 
defines the position of the resonance window and does 
not affect the flipping probability. In contrast, the latter 
is changing with time, and will strongly affect Pwin- 

To account for this effect, we use a mean-field theory, 
treating each spin inside the resonance window equally. 
The interaction is described by adding into the LZ model 
a nonlinear term, Js, where J is the average coupling 
constant between spins and s is the average spin. This 
nonlinear LZ model was first proposed in the context 
of Bose-Einstein condensation iQ. The coupling con- 
stant J is proportional to the average fraction of spins 
in the resonant window and is calculated as J/A = 
Jov/l + 2a/A2i^(-7, - 2Ff, where Jq = di^j) 
is the dipolar field when all the spins point in the same 
direction. 

The probability Pwin is then given by the fiipping prob- 
ability Pniz obtained with this nonlinear model, which 
has been solved numerically and plotted in Fig||. The 
fiipping probability is suppressed for positive J/A and 
increased for negative J/A, compared with the linear LZ 
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probability. Furthermore, we have found that the nonhn- 
ear LZ flipping probability depends on only two parame- 
ters, Pniz = P„;2(a/A^, J/ A), and it has an approximate 
expression O, |l5|. 
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The remaining task is to calculate the distribution of 
local fields. The dipolar field felt by a spin in the lattice, 
consists of two parts: one is the demagnetization field 
from the spins very far away; the other from the neigh- 
boring spins inside a ball Br of radius r ~ (iJu/A)^/'^. 
Since the demagnetization ^dm is contributed by the dis- 
tant spins, it is independent of the lattice structure and 
only depends on the sample shape and the fraction of up- 
spins. Our calculation shows that S^^m = 2CE]j{2F — 1), 
where the constant C is called the shape coefficient and 
can be calculated theoretically On the other hand, 
the dipolar field from the neighboring spins is a distribu- 
tion depending on the lattice structure. With its center 
shifted by the demagnetization field, the overall distribu- 
tion function is 



dk 
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where D{k,F) = Up^^oKl - F)e-''"^'^^^'> + Fe*'=''(^^-)] ||. 
The distribution functions calculated with Eq.p) are 
compared to a Monte-Carlo simulation in Fig. ^there 
is an excellent agreement. 
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FIG. 3: Comparison between our theory and experiments. 
The experimental data on Peg isotopes from Ref.[^, which 
have different tunnel splittings A, collapse on the same curve, 
demonstrating the a/A'^ scaling law. The inset shows the col- 
lapse of the data for Mn4 Ref. |8| , whose A is varied by chang- 
ing the transverse field. The slight deviation in the adiabatic 
regime a/A^ < 0.5 is likely caused by the "hole digging" 
mechanism [ff7|. 



The above discussions indicate that the parameter Jq 
also consists of two parts, Jo = J2b "^(r) + 2CEu- The 



first part only relates to the crystal structure; for the tri- 
clinic Feg, centered tetragonal Mni2, and hexagonal Mn4 
its value is 3.98i?D,1.15ii^D and 12.&2,Ed, respectively. 

Combining Eqs.(l-3), we can integrate the evolution 
equation (|I]) with the initial condition F = 0. The results 
for the tunneling between the two spin states S = ±10 of 
Fcg are plotted in Fig. |[ With the horizontal axis taken 
as a/A^, the experimental data for three different iso- 
topes of Fes collapse onto the curve given by our theory. 
This "collapse" is expected: A and a enter Eq.(|l|) only 
in the combination of a/A^ through Pniz- This remark- 
able scaling law, along with the excellent agreement of 
our theory with the experiments, strongly supports our 
previous argument that nuclear spins do not apprecia- 
bly affect the flipping dynamics of the molecular spins, 
except for renormalizing the tunnel splitting through hy- 
perfine coupling. This scaling law is further confirmed 
by a set of new experimental data on a different system. 
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By changing the transversal mag- 



netic field from OT to 0.085T, the tunnel splitting of Mn4 
is varied almost by an order of magnitude. Nevertheless, 
these data collapse perfectly onto a single curve (inset of 
Fig. |). 

One more important feature in Figs. |l|&^ is the strong 
suppression of the quantum step height, compared with 
the predictions of the LZ model. Two physical mecha- 
nisms are behind the suppression. One is the ferromag- 
netic blocking between spins in the resonance window. 
Due to the sample shape ||l9| and with the shortest lattice 
vector being along the easy axis, the dipolar interaction 
between spins in Fcg is very much ferromagnetic, yielding 
a positive coupling constant J. As seen in Fig. ||, the flip- 
ping probability Pniz is suppressed from P;^ for this case. 
This effect is relatively more significant in the fast sweep- 
ing regime, where there are more spins in the resonance 
window due to the broadened window width and narrow 
dipolar field distribution. Our calculation finds suppres- 
sion of up to 13% due to this mechanism. However, this 
does not account for all the suppression, especially for 
the slow sweeping limit where the resonant window is 
narrower and the distribution function is wider. 

The other mechanism is the shuffling of spins across the 
spectrum of the dipolar distribution D{£^,F). As other 
spins flip, the dipolar field felt by spin i is altered and 
thus gets shuffled to a different part of the spectrum. 
In particular, many spins which are yet to be brought 
into resonance can get shuffled into the swept part of 
the spectrum, losing their chances of flipping. This is 
confirmed by our Monte-Carlo simulation [^Sj, where the 
position of the resonance window is updated after the 
spins in the window are flipped with probability Pwin- 
In Fig. ^, we show how the dipolar distribution function 
changes with the sweeping field in one simulation with 
Pwin — 1- Many spins in the main peak are shuffled 
into the two right peaks. This dominant shuffling to the 
right is related to the largely ferromagnetic character of 
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the dipolar interaction between spins. A careful tracking 
in our simulation shows that about 50% of the spins are 
never brought into resonance and flip zero times, 28% flip 
once, and 12% flip twice. This shuffling mechanism gives 
an intuitive picture of the physics hidden in the evolution 
equation (^. 
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FIG. 4: Change of the internal field distribution in a sweeping 
field. Ed^o ~ l,Awin = 0.1, Pwin — 1. The broadened ver- 
tical line represents the resonance window. The dashed lines 
are calculated from Eq.(^, showing an excellent agreement 
with the solid line. 
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FIG. 5: The dependence of A'/A on the shape coefficient. A 
is the true tunnel splitting; A' is the tunnel splitting measured 
by the LZ method [§. 

With our theory, we can accurately determine the tun- 
nel splitting A of a single molecule from the measured 
quantum step height. The LZ method |j has been used 
to accomplish this, extracting an effective tunnel splitting 
A' from a step height with F/i„ = 1 — exp(— 7rA'^/2a). 
However, the effective splitting A' is not necessarily the 



true splitting A as the LZ model is inadequate to give 
the correct step height. 

Let us consider the fast sweeping limit, where the 
magnetization is very small. In this case, the dipo- 
lar field distribution Eq.(3) is a Lorentzian centered at 
{EB,.d{r) + 2CED){2F-l) with width |^^^(2F)@. 
Then, using the new variables / — F/Piz, x — 'f/EoPiz, 
one can rewrite Eq.lQ) and immediately notice that the 
evolution equation in terms of / and x is independent of 
Piz and E]j. It means that, in the regime a — s- oo, the 
ratio Ff in/Piz = (A'/A)^ tends to a constant depending 
only on C and Jo/Ejy, which describe the geometry of 
the sample: its shape and lattice structure. This asymp- 
totic relation explains why the A' is not necessarily the 
true A of a single molecule. We have calculated the ratio 
A'/A and its dependence on the shape coefficient C for 
three different molecular magnets Fes , Mni2 [jf] , and Mn4 
as shown in Fig. |. With this, one can obtain the real 
tunnel splitting from the corresponding step height given 
the shape of the sample. For the Feg sample used in the 
experiments the shape coefficient C = 1.4||l9|, which 
yields A'/A ~ 0.73. 
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